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We present an update to the NNLL RG-improved QCD prediction of top-antitop production in 
e + e~ annihilation at threshold. It includes for the first time a complete NNLL resummation of 
ultrasoft logarithms, which are dominant at this order and give a sizable correction. The renormal- 
ization scale dependence of the total resonance cross section decreases substantially compared to 
earlier predictions, where the ultrasoft logarithms were included only partially. 
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1. Introduction 

A scan of the resonance in top-antitop production at a future linear collider will allow for a pre- 
cise determination of the mass, the width and the couplings of the top quark This imposes high 
requirements on the precision of the theoretical prediction for the total cross section near threshold, 
which can in principle be met by high order perturbative calculations using nonrelativistic effective 
field theories based on NRQCD [Q]. So far, however the accuracy in the normalization of the total 
1 1 resonance cross section was limited by a relatively large renormalization scale dependence of 
about 6 % in the renormalization group improved (RGI) NNLL prediction ^ (see also Ref. [Q]) 
and about 10 % in a fixed order N 3 LO computation [§J. Although the determination of the top mass, 
which primarily depends on the cm. energy where the cross section rises, might not be affected 
by this uncertainty, it still renders precision measurements of the top total width or the top Yukawa 
coupling impossible. To match the statistical uncertainties, that are expected for these quantities at 
the International Linear Collider (ILC), a theoretical precision of the cross section normalization 
of at least 3% is needed. 

In this talk we shortly review recent progress in the RGI calculations concerned with heavy 
quark-antiquark production at threshold within the vNRQCD effective theory framework [[J and 
discuss its implications for the resonance line-shape of c toi {e + e~ —>ti). 



2. The cross section 



To achieve reliable results for a heavy quark pair production process close to threshold it is 
mandatory to resum so-called "Coulomb singularities", i.e. terms that scale like powers of CC s /v, 
where CC S is the strong coupling and v is the (nonrelativistic) relative velocity of the heavy quarks, 
to all orders in perturbation theory. This task can be performed employing a Schrodinger equation 
within the nonrelativistic effective theory approach. In addition potentially large logarithms of the 
relative velocity v as they typically appear in quantum corrections to the leading order cross section 
can be systematically resummed using a modified renormalization group with the "subtraction ve- 
locity" v being the renormalization scale parameter [|]]. Thus the R-ratio for top-antitop threshold 
production schematically takes the form ^ 

R = -^^=vV ( Y(a,lnvy x (l(LL); a„v(NLL); a. 2 , a s v, v 2 (NNLL) ; ...1 (2.1) 

and we assume v ~ a s due to the Coulombic bound-state-like dynamics of the system in the reso- 
nance region. 

The R-Ratio for the production via a virtual photon or Z boson with the cm. energy y/s has 
contributions from vector and axial-vector currents and can be written as 



where due to the optical theorem 



■ F v (s)R v {s)+F a (s)R a {s) 
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The pref actors F v (s) , F a (s) contain the effects from the y and Z exchange and are given in Ref . |]7|] . 
The standard model currents ft and j" produce the heavy quark pair in a vector and an axial- 
vector state, respectively. Within the effective theory description these currents are replaced by 
their nonrelativistic counterparts and we find to NNLO in the v counting 

471 . 



R v s) 



Im 



s 

4k 



c\ ( v) srf\ (v, m, v) + 2ci ( v) c 2 (v) gf 2 (v, m, v) 



R a (s) = — Im c|(v)^ 3 (v,fn,v) 



(2.4) 



where the =2^ denote effective current correlators and srfi and sfa are suppressed compared to the 
LO correlator 



M(v,m,v) = i £ ■ 
pp' J 



d 4 xe i{ ^- 2m ^ 
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(2.5) 



by a factor of v . The c, ~ 1 are Wilson coefficients of the effective currents, i.e. functions of the 
renormalization scale V and the heavy quark mass m. In Eq. (2^5) we have adopted the vNRQCD 
label notation for the leading 3 Si current 



'LP 



vlo{}<*i)X* 



(2.6) 



where y/ p and Xp annihilate top and antitop quarks with three-momentum p, respectively, and where 
color indices have been suppressed. 

The missing piece to reach NNLL precision in Eq. (2A) is the NNLL running of the S-Wave 
current coefficient ci(v). All other pieces including all relevant electroweak effects 1 (a ew ~ aj) 
are known with sufficient accuracy [?, 10]. In particular, the s$\ correlators can be related to 

Green functions of a Schrodinger operator and incorporate the resummation of Coulomb singular 
terms. 



3. Current renormalization 

Without loss of generality we parameterize the RG evolution of the current coefficient ci(v) 



as 



In 



ci(v) 
ci(l) 



^ NLL (v) + £T LL (v) + r) + . . . , (3.1) 



where the LL contribution vanishes ^ and the NNLL contribution consists of a "mixing" term 
l =m NLL an d a "non-mixing" term <^,^ LL . The latter is generated directly by UV divergent three- 
loop vertex diagrams for the production of a heavy quark-antiquark pair and has been computed 



in Ref. Ql 1| ] using vNRQCD. It involves a non-trivial entanglement of soft and ultrasoft scales at 
subleading order that has at this time only been systematically treated in vNRQCD. The "mixing" 
part, in contrast, is generated by the same two-loop diagrams that are responsibe for the NLL run- 
ning of c\ through subleading NLL corrections to the (four-quark) vertices, namely the potentials, 
in the loops. 



1 Concerning electroweak effects we refer to cross section predictions with moderate cuts on the reconstructed top 
invariant mass []§]. See Ref. for results with more general invariant mass cuts. 
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In other words, we can derive <^ NLL from the known NLL anomalous dimension [JTJ, 13, 14] 



(S 2 = 2) 



V^m[ci(v)] 



NLL 



y c {s) (y) 

16k 2 



?JLM + ( v) + ( V ) + s 2 r W ( y) 



(3.2) 



through the NLL matching and RG evolution corrections to the potential coefficients Y t Jl2| , |i"5| ]. 
Unfortunately the determination of the subleading NLL running of the vNRQCD potentials is not 



yet complete except for the Coulomb potential y} s ^ [|16|, |14j and the spin-dependent potential 



% Q17Q. Recently however we have finished the calculation of the dominant, ultrasoft contri- 



butions to the NLL running of the remaining potential coefficients Y 2 > % and Y^X relevant in 



/t,eff 



Eq. (3.2). In that work [18, 19] we have determined these contributions from vNRQCD diagrams 
with four external heavy quark legs and two ultrasoft gluon loops. Ultrasoft gluons carry energy 
and momentum of Gimv 1 ) and are to be distinguished from soft degrees of freedom with energy 
and momentum of Gimv). The same results were independently obtained also in a different ef- 
fective theory framework called pNRQCD [^(J in Ref. J2"T| ] and allowed us to determine the full 
ultrasoft contribution to the NNLL mixing term <^ NLL in Ref. Jl9|]. 

The analysis of the three-loop (non-mixing) terms in Ref. [ |il| ] showed that the contributions 
involving the exchange of ultrasoft gluons are more than an order of magnitude larger than those 
arising from soft matrix element insertions and in fact similar in size to the NLL contributions. 
The reason is related to the larger size of the ultrasoft coupling a s (mv 2 ) compared to the soft 
coupling a s {mv) and to a rather large coefficient multiplying the ultrasoft contributions. From 
this analysis it is reasonable to assume that the ultrasoft effects, which form a gauge-invariant 
subset, also dominate the mixing contributions. This is also consistent with the small numerical 
effects [^] of the NLL evolution of the spin-dependent coefficient Y s [[n|], which is dominated 
by soft effects and receives ultrasoft contributions only indirectly through mixing. 

The ultrasoft contribution to <^ NLL determined in Ref. [ |l9| ] reads 2 



rNNLL _27tfr A 2 
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(3.3) 



2 The form of Eq. ( |3.3| ) implies the convention that the z and CO — 1/(2 — z) parameters in ^ NLL , as given in Ref. 
are according to Eq. 
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with 

11 4 34 . 20 

Po = —C A --Tn f , p l = —Cl-4C F Tn f -—C A Tn f , (3.4) 

A = -C F (C A + C F )(C A +2C F )—, (3.5) 
S = -C F (C, + C F )(C A + 2C F ) Q(47+ i 6 o ^- 10n ' T , (3.6) 

^(^r^l + ^lnvV 1 . (3.7) 
V a s (m) ) \ 2% ) 

Together with ^,^ LL from Ref. [JTTJ] it gives the dominant contribution to the NNLL running of 
ci(v) (in the MS scheme). The corresponding matching condition ci(l) in Eq. ( fO| ) is obtained 



from a full theory two-loop calculation [£3|, [24p and is given in Ref. [JTTp - 
4. Results 

The numerical analysis of the new result for the NNLL evolution of the current coefficient 
ci (v) carried out in Ref. Jl9| ] shows that the ultrasoft mixing contribution in Eq. (3.3) compensates 



the anomalously large ultrasoft nonmixing contribution to a large extent. That can be observed best 
in the plot of Fig. [I]. The renormalization scale dependence of ci(v) reduces substantially, when 



the new result in Eq. (3.3) is included in the NNLL prediction. This is especially so close to the 
resonance, where typically v ~ v ~ 0.15. 

To illustrate the (small) effect of the known soft contributions at NNLL we have drawn the 
gray area around our new NNLL result (black solid line) in Fig. It shows the variation of the 
result due to multiplying the soft non-mixing terms at NNLL order by a factor between zero and 
two. We believe that this variation represents a good estimate of the uncertainty associated to the 
currently unknown NNLL soft mixing contributions. Compared to the remaining scale dependence 
of the NNLL ultrasoft contributions this uncertainty is negligible. 



Finally we would like to show the effect of the new NNLL result in Eq. ( |3.3| ) on the prediction 
of citot (e + e~ — > 1 1) in the resonance region. In Fig. |awe have plotted the threshold cross section 
at LL, NLL and NNLL order without the new NNLL ultrasoft mixing contribution to the current as 
it was presented in Ref. [^]. In Fig. ^ b the cross section is shown including the new NNLL ultrasoft 
mixing contribution. 3 The decay of the top, actually an electroweak effect, we implemented only 
at LO by shifting the cm. energy yfs — > \fs — iT t . A systematic treatment of electroweak effects 
up to NNLL order can be found in Refs. [§, (26[ |27|, |, [Kj]. 



At any order the velocity renormalization parameter v has been varied between 0.1 and 0.4. 
The reduction of the renormalization scale dependence due to the new NNLL contribution is sub- 
stantial and clearly visible in Fig. || b. However it is conspicuous that the gap between NLL and 
NNLL bands is essentially the same in Fig. |2] a and b and that the remaining NNLL scale variation 
in Fig. § b is comparable to the one of the NLL result. We believe that this is related to an anoma- 
lously small scale dependence of the NLL prediction and might not be interpreted as an instability. 



3 For both predictions the corresponding Schrodinger equation with the NNLL Coulomb potential was solved "ex- 
actly" using numerical methods. Relativistic corrections were consistently treated as perturbations. 
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Figure 1: RG evolution of the 3 5i current coefficient c\ (v) normalized to c\ (1) for top-antitop production 
(m = 175 GeV). The dotted (blue) line represents the full NLL result exp[^ NLL ], the dashed (red) line in- 
cludes in addition the NNLL non-mixing contribution, exp 

[£ NLL + 4iT LL ]- Th e solid (black) line accounts 
for the full NLL and NNLL non-mixing contributions as well as for the new ultrasoft NNLL mixing cor- 
rections, exp[<i; NLL + £,^m LL + C usoftl- The gray area around the black line is generated by varying the 
soft NNLL non-mixing contributions to that curve by a factor between and 2. For the plot we have used 
«i" /=5) (175 GeV) =0.107. 



We postpone the determination of a final number for the remaining theoretical uncertainty of the 
NNLL cross section to an upcoming publication. 



5. Conclusions 

Including the recently obtained ultrasoft NNLL mixing contribution to the running of the lead- 
ing vector current in the prediction for the total production cross section of a top-antitop pair 
near threshold leads to a substantial reduction of the renormalization scale dependence. We ar- 
gue that the uncertainty due to the still unknown NNLL soft mixing contribution to the running 
of the effective current, which is now the only missing piece to a complete NNLL prediction of 
o tot (e + e~ —>■??), can be neglected compared to the size of the now fully known ultrasoft contri- 
butions. These findings will result in a substantial reduction of the perturbative uncertainty of the 
NNLL prediction. A thorough analysis of the respective theoretical errors from all possible sources 
including electroweak corrections and a variation of the vNRQCD matching scale in addition to the 
usual variation of the renormalization scale will be the subject of a future publication. 



6 



NNLL top-antitop production at threshold 



Maximilian Stahlhofen 




342 344 

\/7 (GeV) 



344 346 
\/7 (GeV) 



J4S 



Figure 2: The band labeled "NNLL" represents the prediction of <7 tot (e + e~ —t-ti) in the threshold region 
with (a) and without (b) the NNLL ultrasoft mixing contribution to c\ (v) in addition to all other known 
QCD corrections up to NNLL order. The other bands refer to the complete NLL and LL results, respectively 
and are identical in both panels. All bands were generated by varying the renormalization parameter v 



between 0.1 and 0.4 and using the lS-mass scheme [25 1 with m ls = 172 GeV. Further input parameters 
were T r = 1.46 GeV for the top width and a s (Mz) = 0.1 18 for the strong coupling in the MS scheme. In the 
peak region of the NNLL result in panel b the scale variation is around ± 1 .7%. 
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